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Abstract 

Recently the authors and J.M. Kress presented a special function recurrence relation method 
to prove quantum superintegrability of an integrable 2D system that included explicit construc- 
tions of higher order symmetries and the structure relations for the closed algebra generated by 
these symmetries. We applied the method to 5 families of systems, each depending on a rational 
parameter k, including most notably the caged anisotropic oscillator, the Tremblay, Turbiner 
and Winternitz system and a deformed Kepler-Coulomb system. Here we work out the analogs 
of these constructions for all of the associated classical Hamiltonian systems, as well as for a 
family including the generic potential on the 2-sphere. We do not have a proof in every case 
that the generating symmetries are of lowest possible order, but we believe this to be so via an 
extension of our method. 



1 Introduction 



Recently Tremblay, Turbiner and Winternitz [U [2] studied a family of quantum and classical me- 
chanical systems in two dimensions with classical Hamiltonian 

n I 2 , 1 2 2 2 i a ft ft\ 

n = p r + —p g -ur+- — + , 1 
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We call this the TTW system. They conjectured and gave strong evidence that for k a rational 
number this system is superintegrable, as is the corresponding quantum analog 



H = d 2 r + -d r + -M - 




(2) 



r 2 cos 2 k9 r 2 sin 2 k6 



In papers [3l H] the authors and collaborators proved quantum and classical superintegrability for 
these systems, the first proofs that covered all rational k. Related results are [Hl|6j[7]. This system has 
sparked a great deal of interest because it includes several special cases of definite physical interest, 
it provides many examples of superintegrability where the generating constants of the motion are 
of very high order as polynomials in the momenta or as differential operators in the quantum case, 
and because it gave a hint as to how many other such families of classical and quantum systems 
could be constructed. The previous classical proofs of superintegrability for these systems have 
usually not determined the structure of the symmetry algebras, or even verified closure at finite order 
under commutation. For 2D quantum systems, however, the authors and J. M. Kress have recently 
introduced a method for verifying superintegrability based on recurrence relations for hypergeometric 
functions that enables one to compute the associated structure relations for the symmetry algebra 
with relative ease. In [7] structure relations were obtained for 5 families of superintegrable systems, 
each indexed by a rational parameter k, and including the TTW family. Here we give a simplified 
exposition of some of these results. 

Each of these quantum systems in 2D has an analogous classical counterpart that can be shown to 
be classically superintegrable. Natural questions are 

• What is the structure of the classical symmetry algebra for each of these systems? Does it 
close at finite order? 

• How are the structures of the classical and quantum symmetry algebras, and the symmetries 
themselves, related. 

We will present new results in this paper that provide at least partial answers to these questions for 
each of the systems studied in [7], as well as for a family that contains the 2D system on the 2-sphere 
with generic potential. 

2 Brief resume of definitions 

Suppose we have a Hamiltonian system on a 2D local Riemannian manifold (real or complex) such 
that the Hamilton- Jacobi equation admits additive separation in some orthogonal coordinate system. 
Then the separable coordinates x\,x 2 can always be chosen such that the Hamiltonian takes the form 



u = c 



1 



(Pit +Pl 2 + vifri) +v 2 (x 2 ) 



(3) 
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fi(x 1 )+f 2 (x 2 ) 
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and the phase space function 

is the second order constant of the motion determining the separation. Here, {£2, — where 

2 

is the Poisson bracket of two phase space functions F, Q. A constant of the motion is a nonzero 
function on the phase space whose Poisson bracket with % vanishes. Clearly H and £2 are functionally 
independent constants of the motion, second order in the momenta. We will present strategies for 
determining functions fi, /2, v\, v\ such that there exists a 3rd constant of the motion, polynomial 
in the momenta, and such that the three constants are functionally independent. In this case the 
system is superintegrable. By what is essentially an action angle variable construction |6], one can 
compute the integrals M(x) = \ J y- and N(y) = ~ J ^L. Then C 2 = M — N is a third constant 
of the motion, but usually not a polynomial in the momenta, hence not directly useful in verifying 
superintegrability. In [U E], see also [SI IS], we described a procedure for obtaining a polynomial 
constant from M — N, based on the observation that the integrals 

M= 1 - f-==j± N= 1 - [ / dX2 

can often be expressed in terms of multiples of the inverse hyperbolic sine or cosine (or the ordinary 
inverse sine or cosine), and the hyperbolic sine and cosine satisfy addition formulas. Our method 
works for functions fj,Vj such that M and N possess this property. 

For 2D quantum systems such that the Schrodinger eigenvalue equation = E^fr admits multi- 
plicative separation in an orthogonal coordinate system we have 

H = L 1 = jj-t^jj—M, + 0l 2 + v 1 (x 1 ) + v 2 (x 2 ), (4) 

Jl{Xl) + J2{X 2 ) 

and 

i2 - tsMm <*■ + - uM^ + ^» • 

is the second order operator determining the separation. Here, [L 2 , H] = where [F, G] = FG — GF 
is the commutator of two operators F,G. The separable eigenfunctions ^ = $>(xx)Q(x 2 ) of H are 
then characterized by the equations 

= EV, L 2 m = A\& 

where E is the energy eigenvalue and A is the separation constant. A symmetry operator is a 
nonzero differential operator in xi,x 2 whose commutator with H vanishes. Clearly the second order 
differential operators H and L 2 are algebraically independent symmetries. In [7\ we presented a 
recurrence relation method for determining a 3rd finite order symmetry operator such that the three 
symmetries are algebraically independent, i.e. such that the system is quantum superintegrable. 
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3 The quantum TTW system 



We use the quantum TTW system to illustrate our recurrence relation approach to proving super- 
integrability and determining the structure of the symmetry algebra. For complete details on this 
construction and on its application to other systems, see [7J. The quantum Hamiltonian operator is 
(J2J). The general solution of the eigenvalue problem H^fr = E^fr is 

$ = e -^ T \ k{2n+a+b+l) L k J l 2n+a+b+l \ur 2 )^ (5) 

where we have taken a = k 2 (^ — a 2 ) and (3 = fc 2 (| — b 2 ). The L-functions are associated Laguerre 
and the P-functions are Jacobi, not polynomials in general, [TO]. We consider the functions 

n = e-^ r \ k{2n+a+b+1) L k J 2n+a+b+1 \ujr 2 )P^\x) = Y*{r)X% b {x) 

where x = cos(2k6), X = P or Q and Y = S or T. (Here II is obtained from $ by a gauge 
transformation to remove the angular factors (sin (k6)) a+ 2(cos(k8)) b+ 2.) By this we mean that P%' b (x) 
is a Jacobi polynomial if n is an integer. Otherwise it is given by its hypergeometric expression. If 
X = Q then this denotes the associated second solution of the Jacobi differential equation. Similar 
remarks apply to the choice of Y — S. We have defined this function to be 

and T^(r) to be a second independent solution. The energy eigenvalue is given by 

E = —2lo [2(m + nk) + 1 + (a + b+ l)k] (6) 
and A = k(2n + a + b + 1). The separation equation for 0(0) is 

he = (d 2 + -4— + -4n^)^ d ) = ~ k2 ( 2n + a + b + l) 2 &(6) = -A 2 &(6) 
sm (kO) cos z (kO) 

and L 2 is a symmetry operator for the system. Under the gauge transformation L 2 goes to a symmetry 
that we shall call L 2 and which has the same eigenvalues. We see from the expression for E that in 
order that an energy eigenvalue to be unchanged for different values of m, n we must fix u = m + nk. 
The two transformations 

n — > n + q, m — > m — p 



and 

will each achieve this. 



n — ^ n — q, m — >■ m + p 



Consider the functions X®> (x). If we want to raise or lower the index n we can do so with the 



n 



operators 

(2n + a + b + 2)(1 - x 2 )d x X% b {x) + (n + a + b + l)(-(2n + a + b + 2)x - (a - b))X% b (x) 
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= 2(w + l)(n + a + b + l)^i(x) 

and 

JnX a n \ X ) = 

-(2n + a + 6)(1 - x 2 )9 x X^' 6 (x) - n((2n + a + b)x-(a- b))X% h (x) = 2(n + a)(n + ftJ-X^z), 
[ID] . Similarly, for the functions 3^(i?) = u A l 2 Y^{r) where R = r 2 we can deduce the relations 

(*) = + - j- ^AL4 + 1)} 3tf(i2) = -u^i(^), 



^(5) = { {-A + 1)0* - - + — A(l - A) \ y£(R) = -u{m + l)(m + i4)3&g(fl), 



£ 1 
T + 2iT 

[10J. We now construct the two operators 
and 

When applied to a basis function \I/ n = ^(i^X"' ''l 21 ) f° r fixed w = m + kn, (so m = u — kn) these 
operators raise and lower indices according to 

E+^ n = 2i(-l) p u p (n + l) q (n + a + b+ l) q ^ n+q , (9) 

E_^ n = 2V(-n - a) ? (-n - 6) g (« - A;n + l) p (-« - fc[n + a + 6 + l])p* n - ff , (10) 

where (a) q = (a) (a + 1) - • • (a + q — 1) for nonnegative integer g. It is easy to verify that under 
the transformation n —> —n — a — b — 1 we have E + — >■ E_ and E_ — >• H + . Thus E = E + + S_ as 
a polynomial in n and u is unchanged under this transformation. Therefore it is a polynomial in 
(2n + a + 6 + l) 2 and w. As a consequence of the relation A = —k 2 (2n + a + b + l) 2 , in the expansion 
of S in terms of powers of (2n + a + b + l) 2 and E we can replace (2n + a + b+ l) 2 by L2/k 2 and -E by 
if everywhere they occur, and express E as a pure differential operator, independent of parameters. 
Thus, in the expansion of S in terms of the parameters, a term is replaced by VKif with the W 
operator on the left. Clearly this operator, which we will also call E commutes with H on the formal 
eigenspaces of H . However, in [7] we showed shows that in fact E commutes with if in general, 
thus it is a symmetry operator for if. (Indeed it is shown there that any polynomial operator identity 
established on all formal eigenspaces must hold identically. This is a general fact that holds for all 
applications of the recurrence operator method.) 

We can also easily see that under the transformation n — > —n — a — b — l the operator S + — E_ changes 
sign, hence the operator E = (l/(2n + a + b + 1))(S + — EL) is unchanged under this transformation. 
Again, making the replacements (2n + a + 6 + l) 2 by L2/k 2 and w by —2 — (if + 2uk(a + b+ l))/4u; 
we can express E as a pure differential operator, independent of parameters, and it is a symmetry 
operator for if. Each of these symmetries has a nonzero commutator with L 2 , so each is algebraically 
independent of the set H,L2. This proves that the TTW system is quantum superintegrable for all 
rational k. We set L 3 = S, f 4 = S. 
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Using the explicit relations (Q, ffTU]) for the action of the raising and lowering operators E,± on a basis 
we can obtained very detailed information about the structure of the symmetry algebra generated 
by L 2 ,L 3 ,L 4 . Applying the raising operator to a basis function, followed by the lowering operator, 
we obtain the result 

2-2 + ^ n (11) 
= (-l) p 4 ? w 2p (n + l) q (n + a + l) q {n + b + l) q (n + a + b+ l) q {-u + kn) p (u + k[n + a + b + 1] + l) p * n 

Reversing the order we find 

S+S_^ n (12) 
= (-l) p A q uj 2p (-n) q (-n - a) q (-n - b) q (-n -a- b) q (u - kn + t* -k[n + a + b + l]) p *„ 

Thus the action of the operator + 2 + )2_ on any basis function \l/ n is to multiply it by 

£n + Vn- However, it is easy to check from expressions (jUJ), (fTD] and from f[T2"j) . (fTTj) that under the 
transformations — > —n — a — b—1 we have <H- H + S_ and £ n r/ n . Thus = + 

is an even polynomial operator in (2n + a + 6+l), polynomial in u, and £, n + r] n is an even polynomial 
function in (2n + a + b + 1), polynomial in w. Furthermore, each of and is unchanged 

under the transformation u — > —u — (a + 6+ 1) — 1, hence each is a polynomial of order p in 
[2u + (a + 6+ + l] 2 = E 2 /4u 2 . Due to the multiplicative factor u 2p in each of these expressions we 
conclude that is a symmetry operator whose action on a basis is given by a polynomial operator 
p(+)(F 2 ,L 2 ,w 2 ,a,6). In fact, 

Similarly, the operator S^ - ^ = (H_2 + — H + 2_)/(2n + a+6+l) is an even polynomial in (2n + a + b+l), 
as is (£ n — r] n )/(2n + a + b+ 1). Also it is polynomial in E 2 and w 2 . Thus = P(~\H 2 , L 2 , u 2 , a, b) 
is a symmetry operator which is a polynomial function of all of its variables. 

Now we can compute the structure relations explicitly by evaluating the operators on an eigenfunction 



basis. These relations must then hold everywhere. The results are: 

[L 2 , L 4 ] = -4k 2 qL 3 - 4k 2 q 2 L 4 , (13) 

[L 2 , L 3 ] = 2g{L 2 , L 4 } + 4A; 2 g 2 L 3 + 8k 2 q 3 L^ (14) 

[L 3 , L 4 ] = 2gL 2 - 2P(-\H 2 , L 2 , u 2 , a, 6), (15) 

6fc 2 L 2 + {L 2 , L 4 , L 4 } + 6k 2 q{L 3 , L 4 } + 28k 2 q 2 L 2 4 - 4fc 2 gP ( - ) (# 2 , L 2 , w 2 , a, b) (16) 

-12k 2 p( + \H 2 ,L 2 ,u 2 ,a,b) =0. 

Here, {A, 5} = AL 5 + LM and {A, B,C} is the analogous 6-term symmetrizer of 3 operators. A 
more transparent realization, for R = — 4fc 2 gL 3 — 4/c 2 g 2 L 4 , is 

[L 2 ,L A ] = R, (17) 
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[L 2 , R) = -8k 2 q 2 {L 2 , L 4 } - 16fcV^4, (18) 
[L 4 , i2] = 8k 2 q 2 L 2 4 - 8k 2 qP(-\H 2 , L 2 , u 2 , a, b), (19) 

— + 22k 2 q 2 Ll + {L 2 , L 4 , L 4 } - Ak 2 qP^\H 2 , L 2 , u 2 , a, b) (20) 

-12k 2 P {+ \H 2 , L 2 , u 2 , a, b) = 0. 
From this result we see that the symmetries H, L 2 , L 4 generate a closed symmetry algebra. 

For k = 1 the operator is of order 3, whereas we know that there is a 2nd order symmetry, due 
to the fact that the TTW Hamiltonian-Jacobi equation for k — 1 is multiply separable. Thus we 
have produced only a proper subalgebra of the maximal symmetry algebra. However, for k ^ 1 this 
multiple separation no longer occurs and we think that the generators we have constructed are of 
minimum order. 



4 The classical TTW system 

The TTW Hamiltonian (pQ) admits a second order constant of the motion corresponding to separation 
of variables in polar coordinates, viz 

r 2 a P 

C2=Pe + ^k9 + ^k6- (21) 

To demonstrate superintegrability we need to exhibit a third polynomial constant of the motion that 
is functionally independent of H, C 2 . The method for constructing this symmetry is based on action 
angle variables and the existence of an additively separable coordinate system for the Hamilton- 
Jacobi equation, [U [TT] . We briefly review the construction. In terms of the new variable r = e R the 
Hamiltonian assumes the form 

-2R/J2 i 2 , ,2AR a , P 



n = e- 2 «(p 2 R + p 2 -u J 2 e w + ——- + 



cos 2 (k6) sm 2 (k9) J ' 

Applying the method of [1] to find extra invariants we first need to construct a function M(R,pn) 
which satisfies {M, "H} = e~ 2R , or 

(Au 2 e m + 2He 2R )d PR M + 2p R d R M = 1. 

This equation has a solution 



where 



M = —^=B 

a^/t 2 



. . (2C 2 e- 2R -H) 2^CT 2 e- 2R p R 

smh±5=2 — , coshi)- 



Vn 2 + Au 2 C 2 ' VH 2 + Au 2 C 2 
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and 



, a : (3 



C2 ~ Pd + cos\k9) + sm 2 (key (22) 

and we also have the relation (which we can use to consider M as a function of R alone): 

p 2 R + C 2 -u 2 e m -e 2R 'H = 0. (23) 
We now compute a corresponding function N(9,pg) which satisfies {N, H} = e~ 2R , or 

' " + .,J f ,nJ 9 pB N + 2 P gdgN=l 



K cos 2 (kff) sm\k6) 

where the prime denotes differentiation with respect to 9. This equation has a solution N = ~ 4v ^- fc ^4 

where 

, „ ~P + a - C 2 cos(2k9) , „ y /L 2 'sm(2k9)p e 

sinnA = « — , cosh A - 



y/(C 2 - a - P) 2 ~ 4a/3 y/(C 2 - a - f3) 2 - Aa(3 

According to the action angle variable theory M — N is a constant of the motion and, since it is 
constructed such that {M — N, C 2 } ^ 0, it is functionally independent of C 2 , [T2]. 

From these expressions for M and we see that if k is rational, k — | ( where p, q are relatively 
prime integers) then 

-Aipy^lN - M} = qA + pB 

and any function of qA + pB is a constant of the motion (but not polynomial in the momenta). 
Our original proofs of superintegrability for this system were based on the fact that the hyperbolic 
functions sinh x, cosh x satisfied addition theorems. Here, however, we will make use of an observation 
in [8] that simplifies the analysis and work with the exponential function. We have 

e A = cosh A + sinh A = X/U, e~ A = cosh A - sinh A = X/U, (24) 

e B = cosh B + sinh B = Y/S, e~ B = cosh B - sinh B = Y/S, (25) 

where 



X = V£ 2 sm(2k9)p g + i(-/9 + a - C 2 cos(2k9)), X = V C 2 sm(2k9)p e - z(-/9 + a - C 2 cos(2fc0)), 
Y = 2^/T 2 e- 2R p R + i{2C 2 e~ 2R -H), Y= 2^ 2 e~ 2R p R - %(2L 2 e' 2R - H), 



U = ^(C 2 - a - /3) 2 - 4a P, S = ^H 2 + Au 2 C 2 . 
Now note that e qA+pB and e~( qA+pB ^ are constants of the motion, where 

eqA+P B = {e A )q{e B )P = ^ e - iqA+pB) = {e -A )q{e -B )P = (Xmi_ (26) 

Moreover, the identity e i A +P B e -(i A +P B ) — \ can be expressed as 

X q (X) q Y p (Y) p = U 2q S 2p = P(C 2 , H) (27) 
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where P is a polynomial in £ 2 and H. 



We make some observations. Let a, b, c, d be nonzero complex numbers and consider the binomial 
expansion 

~ ~ ^ ~~ 1 (28) 



£ 2 a + ib) q (^£ 2 c + idf + WC 2 a - ib) q (^£ 2 c - idf 



E 

0<l<q,0<s<p 



{q+p-e~s)/2 



,e+s-\ 



Suppose p + q is odd. Then it is easy to see that the sum (125]) takes the form v / ^2^odd(^2) where 
T dd is a polynomial in C 2 . On the other hand, if p + q is even then the sum (128|) takes the form 
T even (£ 2 ) where T even is a polynomial in C 2 . 



Similarly, consider the binomial expansion 
1 

2i 



E 

0<i<q,0<s<p 



C 2 a + ib) q (^£ 2 c + idf - (^C 2 a - ib) q (^C 2 c - id) p 

£(q+p-e~s)/2 



(29) 



q 



P )b i d s a q - e c p - s ^- 
s I 2i 



Suppose p + q is odd. Then the sum (|29|) takes the form V dd(^2) where V dd is a polynomial in C 2 . 
On the other hand, if p + q is even then the sum (|2"9"|) takes the form v^^Kvent^) where Kven is a 
polynomial in C 2 . 



Case p + q odd: Let 



where 



V ^2 * 



£+ = £~ = (X) q (YY 



(30) 



(31) 



Then we see from (|28|) . ( 129|) that £3,^4 are constants of the motion, polynomial in the momenta. 
Moreover, the identity (|2"T|) becomes 



£+£- = P(£ 2 ,H). 



(32) 



Here £ are constants of the motion, but are not polynomial in the momenta. We will see that they 
are the classical versions of quantum raising and lowering operators for angular momentum. 



It is straightforward to verify the Poisson bracket computations 

jfq-lyp 



{£ 2 ,£ + } = q~i^-{£*,X} = -Mp^C 2 C + 
{£ 2 ,£-} = q { } uq ^ ] {C 2 ,X} = UpyfC 2 C- 



(33) 
(34) 
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From these results we can obtain 



{£ 2 , A} 



Aip^C 2 C + + AipyfC 2 C~ 



-AipC + + AipC~ = ApC 3 



1 



{£ 2 , £3} = -(- WA£ + - Aip^C 2 C~) 
= -4p£ 2 £ 4 . 



To get the remaining structure equations we need some more identities. First of all, note that 

1 



C\ = — [(£+) 2 + 2£+£- + (£-) 2 ] 

A-2 



£ 2 = -[(£+) 2 -2£+£- + (£-) 2 ]. 



Thus, £ 2 + £ 2 £l = 4£+£" = AP(£ 2 ,H), so 



£ 2 = -£ 2 Cl + 4P(£ 2 ,H). 



(35) 
(36) 

(37) 



Further, 



P ( BP BP \ 

{£+, £-} = {£+, — } = (£ + )- 1 {£ + , ^} = (£ + )" 1 ( J^ + > A} + ^{£ + , ^} J 



so 



{£+,£-} =4zp v /£ 2 ~ 



<9P 



(3* 



It remains only to evaluate {£3, £4}. We have 



{£ 3 ,£ 4 } = i{£ + -£-,^= + ^=} 
% V jL 2 v £-2 



- [-(£ 2 )- 3 / 2 {£ + , £ 2 }£ + - (£ 2 )- 3 / 2 {£ + , £ 2 }£- + (£ 2 )- 3 / 2 {£-, £ 2 }£ + + (£ 2 )- 3 / 2 {£-, £ 2 }£- 
z 

+(£ 2 )- 1 / 2 {£ + ,£-} - (£ 2 )- 1 / 2 {£-,£+}] 



Ap 

= ~T 2 

Then, using (|35|) . we conclude that 



[(£+) 2 + 2£ + £- + (£- 



yp 



BP 
d£ 2 



{£ 3 , £4} = -4p£ 2 + 8p 



BP 
dT 2 



(39) 
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To list the structure equations in standard form we define 1Z = 4p£ 3 . Then the structure equations 
take the form 

OP 

{£ 2 , £ 4 } = K, {£ 2 , TZ} = -16p 2 £ 2 £ 4 , {£ 4 , TZ} = !6p 2 Cl - 32p 2 — , (40) 

OL-2 

K 2 = -16p 2 C 2 C 2 4 + 64p 2 P(£ 2 ,H), 

where 

P(£ 2 , H) = ((£ 2 - a - (5) 2 - Aa(5f{H 2 + 4u 2 £ 2 ) p . 
This shows that the symmetry algebra generated by H, £ 2 , £4 closes at finite order. 

Case p + q even: Let 

£ 4 = -^=(C- - £+), £ 3 = £ + + £", (41) 

where £ ± are defined as in fl^Tl) . Then we see from (|28|) . fl2U|) that £ 3 , £ 4 are constants of the motion, 
polynomial in the momenta. Moreover, the identity (|32|) is still applicable, as are (|33|) and ( |34|) . Now 

we obtain 

{£ 2 , £3} = -4^ v / £^£ + + 4ip^/Z~ 2 £- = -ApC 2 d, 

{£ 2 ,£ 4 } = -^={Aip^fr 2 C + Aip^fr 2 C~] = 4pC 3 . 
w£2 

C\ = ~ [(£ + ) 2 - 2£+£- + (£-) 2 ] , (42) 

1-2 

C\ = [(£+) 2 + 2£+£- + (£-) 2 ] . (43) 

we can re derive (j3"7|) again. Using (I3"5|) and (j42]) we re derive (13^1) . Hence the structure equations 
are exactly (1401 again, even though the expressions for £ 3 , £ 4 are different. 

These structure results for classical TTW correspond to the highest order terms in the structure 
equations for quantum TTW, with the correspondence 

*H ~ H, £ 2 ~ L 2 , £3 ~ —kL%, £4 ~ L4. 

Furthermore, modulo multiplication by constants, the raising and lowering operators are associated 
with the classical lowering and raising functions via H± ~ £ T and the quantum recurrences are 
associated with the classical factors via K~ ~ Y, K + ~ Y, J~ ~ X, J + ~ X. Multiplication or 
division of the quantum recurrence operators by 2n + a + b + 1 corresponds to multiplication or 
division of the classical raising and lowering symmetries by y/C^. 



Similarly, since 
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4.1 The classical constant of the motion £ 5 



As for the quantum TTW system with k — 1, our method doesn't give the generator of minimal 

order. In fact, for all rational k, we can extend our symmetry algebra by adding a generator of 

order one less than £ 4 . To see this, consider first the case p + q odd. Then £ 3 = \{C + — C~) is a 

polynomial in £ 2 . Setting £ 2 = in the expressions for X, X, Y, Y we see that the constant term in 
£ 3 is 2(-l)( p -9+i)/2( a _ pyyp Thus 

£ 3 -2(-l)^ +1 >/ 2 (o-/3)^ 

L 5 - -p. (44 J 



is a constant of the motion, polynomial in the momenta. From ()44j) we see that 

£ 2 £ 5 = £ 3 - 2(-l)^ +1 )/ 2 (a - (45) 

so taking the Poisson bracket with £ 2 on both sides of the equation, and using the relation {£ 2 , £3} = 
— 4p£ 2 £zi, we find 

{£ 2 ,£ 5 } = -4p£ 4 . (46) 

Then the identities we have already derived show that £ 2 , £5, Ti generate a closed symmetry algebra 
that properly contains our original symmetry algebra. 

If p + q is even, then £ 3 = £ + + £~. Again £ 3 is a polynomial in £ 2 , but now with constant term 
2(-i)to-*V 2 (a- P) q W. Thus 

_ £ 3 -2(-l)(^)/ 2 ( Q -/3)^ 

^5- -p. I 4 'J 

^2 

is a constant of the motion, polynomial in the momenta. For example, if p — q — 1 then 

C 3 -2(a-/3)H 



£5 



£9 



is a 2nd order constant of the motion. Just as before £ 2 ,£5,"H generate a closed symmetry algebra 
that properly contains our original algebra. 

The construction we have given here is quite general and not at all restricted to the TTW system. 
All that is required is that the constant term in £ 3 be a polynomial in % alone. In all of the examples 
we know, this construction gives the generator of minimum order. However, we don't have a proof 
that the order is minmal for all k. 



5 The classical potential V = k 2 S 2 (x + iy) k 1 /(x — iy) k+1 

Consider the classical Hamiltonian 
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where k — p/q. This system is real in Minkowski space with coordinates X\ — x, yi — —iy. The 
corresponding quantum system has Hamiltonian 



The energy eigenvalue equation is H^f = Ety . In[4j this system was shown to be both classically and 
quantum superintegrable, and in [7] the recurrence relation method was used to obtain the structure 
of the quantum symmetry algebra. Here, we revisit the classical problem and determine the structure 
of the classical symmetry algebra. 

We use the fact that the Hamilton- Jacobi equation is additively separable in polar coordinates 
x = r cos 8, y = r sin 9 where r = e R . As in the previous example it is convenient to pass to variables 
R and 9. Then 

fV 4- n 2 4- k 2 b 2 p 2ike \ 

U = £ 1 = WR + Pe+J de ) } v = k 2 5 2 e 2ike-2 R 
and there is is a constant of the motion 



e 2R 



£ 2 =p 2 9 + k 2 6 2 e 2lke , 

responsible for the variable separation. 

Using the usual prescription for obtaining the extra constant we compute M(R) — \ j N{9) = 

if-- The 

new constant of the motion (not a polynomial) will be M — N. If k is an integer it is 
convenient to consider the solution such that (in terms of the coordinates R, 9) we have 2ipy/£^(M — 
N) = qA + pB as a constant of the motion and 

sinhA = hfi± e -ike CQshA = l JP± e ^e^ 
Ok ok 

iV£~2 -R , „ PR 



sinhi? = — R , cosh B = —=e R . 



m vn 

There is the identity 

e 2R n = p\ + C 2 . 
Using the same reasoning as for the classical TTW case, we have 

e A = cosh A + sinh A = X/U, e~ A = cosh A - sinh A = X/U, (50) 

e B = cosh B + sinh B = Y/S, e~ B = cosh B - sinh B = Y/S, (51) 

where 

X = (ip e + %sJT % )e- ik \ X = (ipe - i^/T 2 )e~ ike , 

Y = {p R + i^T 2 )e- R , Y={p R - i^T 2 )e- R , 

u = k5, s = Vn. 
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Then e qA+pB and e ( qA+ P B ) are constants of the motion, where 

e qA+pB = (e A ne B T = f^, e~^^ = (e^e"*)* = (52) 

The identity e qA+pB e~ ( - qA+pB * > = 1 can be expressed as 

X q (X) q Y p (Y) p = U 2q S 2p = {k5) 2q W. (53) 



Using an analysis similar to (J28J), f l29|) . we see that, independent of p + q being even or odd, we can 
set 

£ 4 = -±={£+ - £"), £ 3 = (£ + + £~), (54) 

V A-2 

where 

£+ = X q Y p , CT = (X) q (7) p . (55) 

Then £ 3 , £ 4 are constants of the motion, polynomial in the momenta. Moreover, the identity ( I53p 
becomes 

£ + £- = (k5) 2q n p . (56) 

It is straightforward to verify the Poisson bracket computations 

Yi- 1 Y P / 

{£ 2 , £+} = g-^-{£ 2 , X} = 2^ v /£ 2 £ + , (57) 

{£ 2 ,£-} = q { } uq ^ ] {C 2 ,X} = -2ip^£ 2 £-. (58) 
From these results we can obtain 

{£ 2 , £ 4 } = - 7 =(2ip^ 2 £ + + 2ip^/T 2 £-) = 2ip£ 3 , 

V ^-2 



{£ 2 , £3} = 2ip\J £ 2 £ + - 2ip^J £ 2 £ = 2ip£ 2 £ 4 . 

Note that 

£\ = _L [(£+) 2 _ 2£+£~ + (£~) 2 ] , £ 2 3 = (£+) 2 + 2£ + £- + (£~) 2 . (59) 
^2 

Thus, 

£ 2 = £ 2 £ 2 + A{k5) 2q W. (60) 

Note that £~ = (k5) 2q W / ' £ + , so 

{£+,£-} = (M) 2 «{£ + , ^} = ^{£ + ,?F} = 
since £ + is a constant of the motion. 
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We find 

{£ 3 ,A} = {£ + + £-,^=-^=} = 

V 1-2 V 1-2 

-(£ 2 )- 3 / 2 {£+, £ 2 }£+ + (£ 2 )- 3/2 {£ + , £ 2 }£~ - (£ 2 )- 3/2 {£-, £ 2 }£ + + (£ 2 )- 3/2 {£-, £ 2 }£~ 

+(£ 2 )- 1 / 2 {£ + , £-} + (£ 2 )~ 1/2 {£~, £+} 

)n, [(£ + ) 2 -2£+£- + (£-) 2 ]. 



£2 

Then, using we conclude that 

{£ 3 ,£ 4 } = 2pz£ 2 . (61) 

To list the structure equations in standard form we define 1Z = —2piC^. Then the structure equations 
take the form 

{£ 2 , £ 4 } = K, {£ 2 , 11} = -Ap 2 C 2 C 4 , {£ 4 , 11} = 4p 2 £ 2 , (62) 

TZ 2 = -Ap 2 C 2 C\ - 16p 2 (kS) 2q W. 

This shows that the symmetry algebra generated by H, £ 2 , £4 closes at finite order. As for the TTW 
case, these classical structure relations correlate with the quantum structure equations in [7] (with 
the sign of the potential reversed) and give the classical version of recurrence relations. 



6 A system on the 2-sphere 

In [7] we considered a quantum Hamiltonian on the two-sphere: 

1 a 

H = d 2 + cotOde + —d% + . 2( , (63) 
sin 6 sin V cos z kip 

where k = | and p, q are relatively prime positive integers, showed it was superintegrable, and 
determined the closure of its symmetry algebra via the recurrence relation approach. Now we treat 
the classical analog 

K = V 2 e + -K- ft vl+ . 2fl Q 2 , • (64) 
sin T sin 6 cos z kip 

In this case the variables are not in the standard form for our classical construction and we need 
to consider the separation variables x±,x 2 as </?, if? where x% = ip, x 2 — ip, fi = , f 2 = 1/ cosh 2 ^, 
V\ = a I cos 2 kip. Then 



a 



U = cosh 2 ^(pl +pl + — rT r) (65) 
v r cos kip 



and 



a 



C * = ti + ^nri ( 66 ) 

* COS z 
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so H = cosh 2 ipijpQ + C 2 ). The variables ip,9 and their corresponding momenta are related by 
cosh-?/; = 1/ sin^, sinh-0 = cot 9 and p^ = — sin 9 p g . Then 



1 f dip 1 



2J P v 2 



so 



and 



M 



-%A 



£2 tit- ^ 2 

^ cos^ fe(/5 



, sinhA=^ SinA:V9 , coshA= C -^m 



dip 



n -£2 



cosh 2 ip 



\/C 2 — a 



N 



iB . — iky/C^ sinh-0 —ik\fll 2 cot9 . 
, smhiJ = — = — , coshi? 



2v/£ 



^U-C 2 



^U-C 2 



\/C 2 — a ' 

cosh V> _ -pg 



From these expressions for M and N we see that 

1ip\/T 2 \M -N] = qA + pB 

is a constant of the motion. As usual, 

e A = cosh A + sinh A = X/U, e~ A = cosh A - sinh A = X/U, 
e B = cosh B + sinh 5 = Y/S, e~ B = cosh 5 - sinh B = Y/S, 



(67) 
(68) 



where 



X = cos ktp p v + i \/C 2 sin ktp, X = cos ktp p v — i \/C 2 sin kip, 
Y = — pq — ik\l~L 2 cot 9, Y = — p e + ik\pi 2 cot 9, 
U = ^C 2 - a, S = ^/U - C 2 . 



Again, e qA+pB and e ( <?j4+pB ) are constants of the motion, where 



^qA+pB 



(e A ) q (e B ) P = e -{lA+pB) _ ( e ~A\q( e -B^p 



(X)i(Y)p 



and the identity e <i A +P B e -(<i A +P B ) — \ can be expressed as 

X q (X) q Y p (Y) p = U 2q S 2p = P{C 2 ,U) = (C 2 - a) q {U - C 2 ) p 
where P is a polynomial in C 2 and H. 



U q S p 



(69) 



(70) 



As in the example immediately preceding, we need consider only one case, independent of p + q even 
or odd: 

£ 4 = -^={C + - £-), £ 3 = C + + C~, (71) 



where 



£+ = X q Y p , cr = {X) q {Yf. 



(72) 
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Then £ 3 , £ 4 are constants of the motion, polynomial in the momenta, and 

c + c- = p(c 2 , u) = (£ 2 - a y(u - c 2 y. 

Now, 

Xi- l Y p , — (X) q ~ 1 (Y) p , — 

{£ 2 ,£+} = q— — {£ 2 ,X} = -2ip^T 2 C\ {C 2 ,C-} = q [ ) {C 2 ,X} = 2ip^C 2 £ 



UiSp 



UiSp 



so 



Since 



-2ip^C 2 C + - 2ipy/£ 2 £-) = -2ipC 3 , 



{£2, £4} 



{£ 2 ,£ 3 } = {-2ip^C 2 C + + 2ip^JC 2 C-) = -2ip£ 2 £ 4 . 



£ 2 = — [(£+) 2 - 2C+C- + (£-) 2 ] , £ 2 = (£+) 2 + 2£+£~ + (£~) 2 , 



we have £ 2 - C 2 C\ = 4£+£- = 4P(£ 2 , U), so 

£ 2 = £ 2 £ 2 + 4P(£ 2 ,?/). 



Further, 



P f BP BP \ 

{£+, £-} = {£+, — } = P} = CC+)- 1 [g^i^, A} + m {C\ n}j , 



so 



{£+,£-} = 2%p^/T 2 



dP 



Lastly, 



{£ 3 , A} = {£+ + £-, 4^ " 4^} 



-(£ 2 )- 3 / 2 {£+, £ 2 }£+ + (£ 2 )- 3 / 2 {£ + , £ 2 }£~ - (£ 2 )- 3/2 {£~, £ 2 }£ + + (£ 2 )- 3/2 {£~, £ 2 }£' 

-(£ 2 ) _1/2 {£ + , £"} + (£ 2 ;r 1/2 {£~,£ + } 



= 2ip£l — Aip 



dP 
dTo 



For standard form we set 1Z = —2ip£ 3 . Then the structure equations take the form 



{£ 2 , £4} = n, {£ 2 , n} = -a p 2 c 2 c 4 , {£ 4 , n} = -A P 2 ci + 8 P 



dP 



H l = -4p 2 C 2 C 2 4 - l6p z P{£ 2 ,U), 



where 



P(C 2l U) = (C 2 -af(H-C 2 y. 
Thus the symmetry algebra generated by H, £ 2 , £4 closes at finite order. 
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7 The caged anisotropic oscillator revisited 



The quantum caged anisotropic oscillator is defined by the Hamiltonian operator 

H = d 2 x + d 2 y - uj\x 2 - u 2 2 y 2 + 4 + 4' ( 8 °) 

x y 

where U\ = pu and u> 2 = qoo and p, q are positive integers that we assume are relatively prime. The 
classical version is 

IT 2,2 22 2 2 , a , @ /oi\ 

H = P x + P y -u 1 x -u 2 y + — + —, (81) 

x y 

Marquette [TH [To] , using a ladder operator method, was the first to work out the quantum structure 
relations for this system. We treated it again, from the recurrence relation point of view, in [7] and 
will revisit it here to demonstrate the correspondence between the quantum recurrence relations and 
our classical construction. We do not claim novelty for the final results for this very simple but very 
important system, only the approach. This system separates in the Cartesian coordinates x, y and 
either of the classical constants of the motion 

r 2 2 2 i a r 2 2 2 i P 

£i=Px~ u x x + — , £ 2 =p y - u 2 y + — , 

x y 

characterizes the separation, where H = Ci + C 2 . 

Constructing M(x,p x ) which satisfies M — ~ J dx/p x , we find a solution 

M = — , sinh A = ^ , cosh A - ^ + A 



W ^Cj + Au 2 a y/Ll + 4u>& 

Similarly, N(y,p y ) can be found as 

N = sinh B = , cosh B - ^ + ^ 



W y/L 2 2 + 4u*p' v^f+4^p' 

Thus, 4pqu)(M — N) = qA + pB is a constant of the motion, as are 

e« A+pB = (e A ne B r = f^, = (e^e^ = (82 ) 

and the identity e i A +P B e -(i A +P B ) — \ can be expressed as 

X q (X) q Y p (Y) p = U 2q S 2p = P{C Xi C 2 ) = {C\ + 4u; 2 a) 9 (£ 2 + Au 2 a) p (83) 
where P is a polynomial in C\ and In particular, 

X = 2x 2 w 2 + L\ + 2a; 1 xp a ;, X = 2x 2 w 2 + £i — 2w!a;p x , 
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Y = 2y 2 u 2 + £ 2 - 2w2VP v , Y = 2 Z/ 2 ^2 + £2 + 2u 2 yp y , 
U = ^Lj + 4oo 2 a, S = ^L 2 + 4oo 2 /3. 

As usual, we set 

£+ = X q Y p , £- = (X) q (Y) p . (84) 

In this degenerate case, £ + , £~ are each constants of the motion, polynomial in the momenta, We 
set 

£ 3 = £+ + £-, £ 4 = £+-£-. 

Now 

{d,C + } = qX^Y^X} = -4g Wl /: + , {&,£-} = Aq^C , 
{£ 2 , £+} = pX q Y^ l {£ 2 , Y} = 4ptu 2 £+, {£ 2 , £'} = -Apu 2 £^ 
and, since C + C~ = P, 

{£+,£-} = {£\^} = (£ + )- 1 {£\P} 
= (£ + r l ( A} + SH£ + , £2}) = 4g Wl |^ - A P u 2 9P 



d£^ ' ±J d£ 2 l ' "J " A ad " ^<9£ 2 
Thus, 

{£i,£ 3 } = -4gwi£ 4 , {£i,£ 4 } = -4gwi£ 3 , {£ 2 ,£ 3 } = Apu 2 £ 4 , {£2, £4} = Apu 2 £ 3 , 

dP dP 
{£1, £ 2 } = 0, {£ 3 , £4} = 8gwi— - - 8pw 2 — -. 

a£i a£ 2 

£\-£\ = 4P, P(A, £ 2 ) = (£ 2 + 4 Wl 2 a)' ? (£ 2 + 4 W 2 a) p . 

8 Extension of the generic potential on the 2-sphere 

A general Hamiltonian on the 2-sphere, separable in spherical coordinates can be taken in the stan- 
dard form 

U = cosh 2 ^){p\ + pi + v^ip) + v 2 (ip)), x = <p, y = 4>, fi = 0, f 2 = — - . 

Embedded in complex Euclidean 3-space with Cartesian coordinates, such 2-sphere systems can be 
written as 

H = j* + j2 + ji + v(s), 

where 

Js = SiP S2 - s 2 p Sl , Ji = s 2 p S3 - s 3 p S2 , J 2 = s 3 p si - Sip S3 , s\ + s\ + s\ = 1. 
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where s± = sin#cos</9, s 2 = sin^sin^, S3 = cos# are the usual spherical coordinates. The variables 
ij},9 and their corresponding momenta are related by cosh?/' = 1/ sin#, sinh^ = cot# and p^ = 
— sin# pq. Now we consider the superintegrable Hamiltonian system 



H = cosh 2 i,\p\ +pl + + + — ' ] (85) 

r cos z kip sm kip smh ip 

£ 2 a ft 

2 ^ v cos 2 kip sin 2 fc<^ 

Here k = p/q where p, q are relatively prime integers. As shown in [7], the functions that determine 
the extra constant of the motion are 

MfapJ = — i -i=A sinh^ = i(f2Cos(2k<p)-a + P) ^ v^sin(2^K 



±ksfT 2 ' v/(£ 2 -«-/3) 2 -4a/3' " ^(£3 - a - /3) 2 - 4a/?' 

AT(V^) = -Uj, sinhi? = y2Cosh(2V>)+ 7 -^) ? ^ _ ^sinh(2^ 



4v/A ' v/CH-£ 2 -7) 2 -4£ 27 ' ' ' v/CH-£ 2 - 7 ) 2 -4£ 27 ' 

(Here, we have corrected some misprints in the expressions listed there.) If k = 1 this is the generic 
potential on the 2-sphere, a 2nd order superintegrable system. 

Carrying out our usual construction, we have 

e A = cosh A + sinh A = X/U, e" A = cosh A - sinh A = X/U, (86) 
e B = cosh B + sinh 5 = Y/S, e~ B = cosh 5 - sinh B = Y/S, (87) 

where 

X = ^/C 2 sin(2k(p)p (f + i((3 — a + C 2 cos(2fc<p)), X = yX^sin^A;^)^ — z(/3 — a + C 2 cos(2kip)), 

Y = v / Z^sinh(2V')^ + i( 1 -H + C 2 cosh(2^)), Y = v / ^sinh(2^)^ - 1(7 - H + £ 2 cosh(2^)), 
C/ = v 7 ^ - « - /3) 2 - 4a/3, 5 = v 7 ^ - £ 2 - 7 ) 2 - 4 7 £ 2 . 
Then _ _ 

e« A+pB = (e A Y(e B r = = (e^e"*)* = (88) 

and we have the identity 

X^Z^Y^Ff = f/ 2 ^ = P(£ 2 , H) = [(£2 - a - /3) 2 - 4a/?] 9 [(ft - C 2 - 7 ) 2 - 4 7 £ 2 f (89) 

where P is a polynomial in C 2 and "H. 

Case p + q odd: Let 

£ 4 = ^=(£ + + £-), A = -(£ + - £-), (90) 
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where 

£+ = X q Y p , CT = (X) q (Y) p . (91) 
Then £ 3 , £ 4 are constants of the motion, polynomial in the momenta and we have 

C + C- = P(C 2 ,H). (92) 

We find 

vq-iyp , 

{£ 2 , £+} = g-^-{£ 2 , X} = 4^v^£ + , (93) 

{£ 2 ,£-} = g Uq g P i^,X} = -Mp^C 2 C-, (94) 

so 

{£ 2 , £ 4 } = -4p£ 3 , {£2, £3} = 4p£ 2 £ 4 . 



Using f l35|) . (j36j) we obtain 

Cl = -C2Cl + 4P{C2,H). (95) 

Further, 

P /r)P <9P \ 

{£+, £-} = {£+, — } = (£ + )- 1 {^ + , P} = (C + )~ l (^H^* ^ + ^{£ + , n}j , 

so 

/ — dP 

{£+,£-} = -4ip^/Zr 2 —, (96) 

and 

dP 

{£ 3 ,jr 4 } = 4pjr,l-8p—. (97) 

To list the structure equations in standard form we define 1Z = —ApC^. Then the structure equations 
take the form 

dP 

{£ 2 , £4} = n, {£ 2 , 11} = -16p 2 £ 2 £ 4 , {£ 4 , TZ} = lQp 2 C\ - 32p 2 -— , (98) 

OL 2 

TZ 2 = -16p 2 £ 2 £ 2 + 64p 2 P(£ 2 ,U), 

where 

P(C 2 , U) = ((£ 2 - a - f3) 2 - Aaf3) q ({H - £ 2 - 7 ) 2 - 4 7 £ 2 )*>. 
This shows that the symmetry algebra generated by "H, £ 2 , £4 closes at finite order. 

Case p + q even: Let 

£ 4 = r4=(£- - £3 = £ + + £", (99) 
w£2 

where £ ± are defined as in ( 19T|) . Then £3, £4 are polynomial constants of the motion. It can be 
verified that the structure equations are exactly the same as for p + q odd. 
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Just as for the TTW system, in the case k — 1 our method gives a 3rd order generator, whereas 
we know that, due to the multiseparability of the Hamilton- Jacobi equation, there is a 2nd order 
generator £ 5 independent of L^^l-i- The method for finding this symmetry £ 5 for all rational k and 
showing that it determines a larger symmetry algebra that closes is exactly analogous to that in 
Subsection 14.11 The message is that the symmetries C determine everything. 

9 Discussion and conclusions 

We have found the classical analog of the recurrence relation approach for verifying quantum superin- 
tegrability for 2D Hamiltonian systems that admit a separation of variables, and have demonstrated 
that this approach allows us give explicit expressions for the symmetries and to compute structure 
equations for them. We have not proved that this method always yields the symmetry generators of 
lowest order, although that is the case for all of the examples we have checked. A major surprise for 
us was the simplicity of the structure equations. Indeed all of the structures in our various examples 
are essentially the same, except for the function P(£ 2) %) and the polynomial in % that occurs in 
the expression for £5. It is clear that this method can be generalized to Hamiltonian systems in > 2 
dimensions that admit separation in subgroup type coordinates [3]. It is not clear what to do if the 
Hamiltonian system does not admit orthogonal variable separation. 
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